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Abstract 

At high temperature the infrared modes of a weakly coupled quan- 
tum field theory can be treated nonperturbatively in real time using 
the classical field approximation. We use this to introduce a nonper- 
turbative approach to the calculation of finite-temperature spectral 
functions, employing the classical KMS condition in real time. The 
method is illustrated for the one-particle spectral function in a scalar 
field theory in 2 + 1 dimensions. The result is compared with re- 
summed two-loop perturbation theory and both the plasmon mass 
and width are found to agree with the analytical prediction. 
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1 Introduction 



Finite-temperature field theory has received considerable attention during recent 
years (see for a comprehensive textbook). An important motivation is the 
physics of the quark-gluon plasma, currently under investigation at RHIC, as are 
baryogenesis and reheating after inflation in the early universe. Thermal field 
theory also provides a necessary reference point for the more complicated case of 
nonequilibrium quantum fields. 

In thermal equilibrium a prominent role is played by spectral functions since 
other correlators can be recovered from it via the Kubo-Martin-Schwinger (KMS) 
periodicity condition |2j . Thermal field theory problems can therefore be reduced 
to a calculation of the appropriate spectral function. In particular the one-particle 
spectral function contains information on the quasiparticle structure of the the- 
ory, needed to describe transport properties of hot matter in the quark-gluon 
plasma with the help of a Boltzmann equation ||. For the calculation of trans- 
port coefficients, such as the shear viscosity, the presence of a medium-dependent 
finite width in the one-particle spectral function is crucial Q. Resummed pertur- 
bative descriptions of the equation of state of the hot QCD plasma may require 
a consistent inclusion of nontrivial quasiparticle spectral functions |J. 

In spite of the apparent importance of spectral functions, nonperturbative 
computational schemes are rather scarce. A first-principle approach is offered by 
lattice field theory. However, the necessity to use a euclidean formulation hinders 
access to dynamical quantities such as spectral functions and other real-time cor- 
relators. Experience in the recovery of mesonic spectral functions in QCD from 
euclidean-time correlators has been gained in the last few years using the Maxi- 
mum Entropy Method (MEM) ||. At high temperature this approach becomes 
especially difficult due to the compactness of the euclidean-time direction. A 
formulation directly in real time avoids these problems. Unfortunately a fully 
nonperturbative approach to real-time quantum field correlators is still lacking.^] 

At high temperature and weak coupling a nonperturbative approach to real- 
time quantities is provided by the classical approximation, originally proposed in 
the context of high-temperature sphaleron transitions and electroweak baryogen- 
esis |§ . Indeed, at high enough temperature the infrared sector of a thermal field 
theory behaves classically as can be guessed from the Bose-Einstein distribution 
function at low spatial momenta: 

1 T 

n(uJ P ) = 7T 77^ t -> t — = n cl (u p ), hup < T, (1) 

exp(n-a;p/T) — 1 nujp 

with ujp = y/p 2 + m 2 and T the temperature (we take H — 1 from now on). As 
is well known, a proper definition of classical thermal field theory requires an 

1 Recent progress in nonequilibrium quantum field dynamics using the 2PI effective action, 
including a calculation of the out-of-equilibrium spectral function for a scalar field in 1 + 1 
dimensions, can be found in jTj. 
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inherent ultraviolet cutoff, provided for instance by a spatial lattice, to regulate 
the Ray leigh- Jeans divergence. The importance of the interplay between the 
ultraviolet lattice modes and the physical infrared modes has been realized first 
in Ref. ||. Much progress in the understanding of the classical approximation 
and quantum and classical thermal field theory has been made subsequently, both 
numerically ]10| and analytically ||TT| , [12| , culminating in Bodeker's effective 
theory for hot infrared nonabelian field dynamics [Q. A recent review discussing 
various aspects of the classical approximation can be found in Ref. JT5|| . 

In this paper we introduce a nonperturbative approach to the calculation of 
spectral functions using the classical field approximation (Sec. 2). We demon- 
strate the method with a calculation of the one-particle spectral function in a 
scalar field theory in 2 + 1 dimensions in Sec. 3. In Sec. 4 we calculate the re- 
summed perturbative spectral function and contrast it with the nonperturbative 
numerical result. Our findings are summarized in Sec. 5. For a discussion of the 
classical analogue of thermal field theory for a weakly coupled scalar field we refer 
to Ref. ra. 



2 Classical approximation 

We consider an arbitrary bosonic operator O and define the spectral function as 
i times the expectation value of the commutator 

p(x- y )=i([0(x),Oi(y))J). (2) 
The brackets denote expectation values at finite temperature T, 

(O) = ^Tr e- H / T 0, Z = Tr e~ H / T , (3) 

where the trace is taken over the Hilbert space. Operators 0(x) = 0(t, x) are 
time dependent with the time evolution determined by the hamiltonian H, 

<3(t,x) = e iHt O{0,x)e-' lHt . (4) 

The spectral function obeys p^(x) = — p(— x) and is, in our convention, real for 
an hermitian operator: 

O 1 " = O -> p\x) = p(x). (5) 

In equilibrium two-point functions depend on the relative coordinates only and 
it is convenient to go to momentum space, 

p{p) = J d 4 xe- ip - x p(x), (6) 
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where p ■ x = —p°x° + p • x, p° = uj = E and x° = t. We find that in mo- 
mentum space p{p) = ipi m (p) is purely imaginary. The imaginary part obeys 
P°Pim(p°,p) > and p im (-p) = -Pimip). 

For a straightforward discussion of the classical approximation it is convenient 
[ p~3|j to introduce also the Keldysh or statistical two-point function fT6|j , 

F(x-y)= 1 -([0(x),0\y)} + ), (7) 

obeying F\x) = F(—x), F^(p) = F(p). In terms of the usual Wightman func- 
tions ffl, 

G>(x -y) = (O(x)Ot(y)), G<(x -y) = (0\y)0(x)), (8) 
the spectral and statistical two-point functions read 

p(x) = i [G>(x) - G<(x)) , F(x) = \[G>{x) + G<(x)}. (9) 

In equilibrium the importance of the spectral function is manifest since all 
two-point functions introduced above can be expressed in it, due to the KMS 
condition. We find in particular |L3 



F(p) 



«(P°) + 1} P(P), ndf)= 1 ^ (10) 



where n(p°) is the Bose-Einstein distribution function. This relation is exact. The 
question is whether the spectral function can be computed nonperturbatively. 

The classical approximation allows access to nonperturbative correlation func- 
tions in real time, both in and out of thermal equilibrium. In a classical theory 
operators commute and the basic classical equilibrium correlation function in a 
bosonic (scalar) theory is given by 

S(x-y) = (0(x)0\y)) cl =^- [ DttD^ e- H / T 0{x)0\y), (11) 

with the classical partition function Z c \ = J D7rD<f) exp(—H/T) and 0(x) = 
0[<f)(x),7i(x)]. This correlator is the classical equivalent of the Keldysh two- 
point function ([?[). The functional integral is over classical phase-space at some 
(arbitrary) initial time, 

J DnD(j) = J I]d7r(x)d0(x), (12) 

weighted with the Boltzmann weight, providing initial conditions 0(0, x) = 0(x) 
and 7r(0, x) = 7r(x). The subsequent time evolution is determined from Hamil- 
ton's equations of motion for 4>(x) and tt(x). The definitions are given for a scalar 
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field theory, but they can easily be carried over to (non)abelian gauge theories 
plj. The most convenient formulation employs the temporal gauge, with the 
Gauss constraint imposed on the initial conditions. It is subsequently preserved 
by the classical equations of motion. 

The classical spectral function is obtained by replacing — i times the commu- 
tator in Eq. (JJ]) with the classical Poisson brackets, 



Pc i(x-y) = -({0(x),0\y)}) 



cl- 



defined with respect to the initial fields 



{A(x),B(y)} 



d 3 : 



5A{x) SB(y) 8A(x) 5B{y) 



5(f)(z) §n(z) 8ir(z) S(j)(z) 



(13) 



(14) 



Due to the formal correspondence between commutators and Poisson brackets 
the quantum and classical spectral function obey the same basic properties. 

At first sight a calculation of the classical spectral function from the defini- 
tion in terms of the Poisson bracket appears rather hard. Fortunately, in thermal 
equilibrium we may use the KMS condition to simplify the procedure. The clas- 
sical KMS condition is based on the same principle as the usual KMS condition 
in a quantum theory: the thermal Boltzmann weight and the time evolution are 
controlled by the same hamiltonian H . An easy way to find the classical KMS 
condition is to consider the high-temperature (or h — > 0) limit of the quantum 
KMS condition. The classical equivalent of Eq. (jlOf) reads (compare Eq. ([!])) 



T 

S ^ = _i ^d(p). 



(15) 



One may also derive this relation directly in the classical theory without refer- 
ence to the quantum case 



17, [L3|. This leads to the classical KMS condition 



formulated in real space, 



p d (t,x) = -—dtS(t,x). 



(16) 



This relation will form the basis of the remainder of this paper. It allows us 
to calculate the spectral function in real time from a relatively easily accessible 
unequal-time correlation function. 

3 Scalar field on the lattice 



As an example we discuss the simple case of a real scalar field in 2+1 dimensions 
with the hamiltonian 



H 



d x 



i 7r 2 + I (V0)2+ i 2 2 + A 4 

2 2 v ; 2 ^ 4! 



(17) 
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Figure 1: Classical spectral function p c \(t, 0) at zero momentum in real time. The 
inset shows the absolute value on a log scale. The temperature is T/m = 7.2. 
Without loss of generality we use X/m = 1 throughout. 



We focus on the one-particle spectral function and take 0(x) = 0\x) = (f>(x) 
Using Eqs. (|16j) and (|TT| ) we find directly 

Pd(t,x) = --(7r(*,x)0(O,O)) d . (18 



The right-hand-side of Eq. fll8|) can be computed numerically in a straightforward 
manner, as follows. The classical field theory is defined on a spatial lattice with 
N x N sites and lattice spacing a and we use periodic boundary conditions. To 
solve the dynamics we use a leapfrog discretization with time step ao < a. The 
canonical momenta 7i(t + |ao, x) are defined at intermediate time steps, which 
suggests to use a symmetrized definition of the spectral function on the lattice 

Pci,iat(*,x) = -^(?r(t + ^o,*)^ [0(0,0) + <P{a ,0)]\ (19) 
To generate thermal initial conditions we use the Kramers equation algorithm, 



a variant of the hybrid Monte Carlo method |p!8| . The evolution in real time is 
calculated using classical equations of motion. In a simulation we switch there- 
fore between noisy evolution to create independent thermal configurations and 
hamiltonian evolution to calculate observables. The results presented below are 
obtained using 2000 independently thermalized initial configurations for each 
temperature. The mass scale m is used as the dimensionful scale and the results 
presented are obtained with N = 128, am = 0.2 and a /a = 0.1 (note that the 
finite time step affects the equal-time canonical relation d t p c i(t, x)| t=0 = 5(x)). 
In a classical theory the coupling constant A can be scaled out of the equations 
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Figure 2: Spectral function p c \(uj, 0), obtained from the real-time result of Fig. [I] 
by a sine-transform, versus the frequency uj. The inset shows a magnification with 
the data points indicated. The dotted line is a fit to a Breit-Wigner function. 



of motion and the remaining dimensionless combination is XT/m 2 (recall that A 
has a dimension of mass in 2 + 1 dimensions). Without loss of generality we take 
therefore X/m = 1 throughout. Larger T corresponds then to a larger effective 
interaction strength. In the simulations the temperature is determined from the 
average kinetic energy T = a 2 (n 2 (t, x)) c i and the temperatures we consider are 
such that aT = Oil). 

In Fig. 1 we present the classical spectral function at zero spatial momentum, 
obtained from a volume average of Eq. (|19D , in real time at a temperature T/m = 
7.2. We see oscillating approximately exponentially damped behaviour. The 
spectral function in frequency space can be obtained from a sine-transform 

p t max 

p c i(a;,p) = 2/ dt p c i M (t,p)sinujt, (20) 
Jo 

where we used the antisymmetry of the spectral function under time reflection 
and absorbed the i of the previous section directly in the definition. The result 
is shown in Fig. 2. The spectral function consists of a single peak with the 
narrow width. The inset shows that the peak is well described by a Breit-Wigner 
spectral function (see Sec. 4). We have tried to find other contributions at larger 
frequency, but these could not be seen in the numerical data. In the numerical 
simulation the late-time regime becomes more and more difficult to establish, 
resulting typically in wiggly nondamped behaviour. This limits the time interval 
that can be used in the sine-transform to a maximal time t max and constrains 
the resolution in frequency-space to Au; = 7r/i max . For the result in Fig. 2 the 
resolution is Au/m = rr/300 ~ 0.01, as can be seen from the inset. 
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4 Perturbative expectation 



The spectral function in the quantum theory can be expressed in terms of the 
retarded self energy = Re £r + 2 Im E# as JO^ ^(J 

/ n_ -2ImE R (^,p) 

[u; 2 - p 2 - m 2 - Re p)J + [Im £r(u;, p)J 

The theory exhibits a quasiparticle structure, with the quasiparticle often referred 
to as the plasmon, in the limit that the rate T(u;,p) = -ImSfi(w,p)/u is much 
smaller than (p 2 + m 2 + ReE^) 1 / 2 . In this case the spectral function is well 
approximated with a Breit-Wigner function, which at zero momentum reads 

2uT 

Pbw(w ' 0) " (c , 2 _M 2 ) 2 + u; 2 r 2 - (22) 



Here M is the plasmon mass and V its width (at zero momentum) 

ImS fl (M, 0) 
M 



(23) 



In this limit contributions from multiparticle states beyond the three-particle 
threshold are tiny. 

A perturbative calculation of the retarded self energy is standard in thermal 
field theory (see e.g. for a clear discussion in 3 + 1 dimensions). For the 
(2 + l)-dimensional case we consider here we find the following. At one-loop 
order the tadpole diagram 

^ = ±[p^A±A (24 ) 

R 2 J (2tt) 2 uj p k 1 

contributes to the mass shift only. Resummation of the tadpole diagram in the 
limit of high temperature and weak coupling results in a gap equation for the 
resummed mass parameter M: 

XT T 

M 2 = — log — (one-loop resummed) , (25) 
Air M 

where the zero-temperature mass is neglected. In 2 + 1 dimensions the one-loop 
mass is sensitive to both the ultraviolet momentum scale (cutoff by T) and the 
infrared momentum scale (cutoff by M). A finite width in the spectral function 
arises at two-loop order from the imaginary part of the setting-sun diagram. We 
focus here on on-shell 2^2 scattering for which the contribution reads 

A 2 r 

ImE^Wp.p) = — —J d$i 23 (p) 2tt5(^ p + co k - tu q - u T ) 

x {n(uj k )[l + n(uq)][l + n{u T )\ - [1 + n(^ k )]n(^ q )n(u; r )} , 



S 




Figure 3: Spectral functions p c i(u,0) for various temperatures T. Fits to a 
Breit-Wigner function are shown with dotted lines. The height, position and 
width of the peak are approximately related as p c \(M, 0) ~ 2 /MY. 

with 

d< W P ) = (9 d l (9 d 'l (f) % (2 7 r) 2 ^(p - k - q - r). (26) 
(27r)^2o;k {2k) lu)^ (27c) z 2uj r 

The on-shell dispersion relations contain the one-loop resummed mass parameter, 
uj p = y/p' 2 + M 2 . It is easy to check that the momentum integrals are dominated 
by the infrared modes. These soft modes are classical. The leading contribution 
at high temperature and weak coupling can therefore be obtained by replacing 
n{uj) — > T /uj = n c \(uj). At zero momentum p = the integrals can be performed 
analytically and after a straightforward calculation we find the plasmon width in 
2 + 1 dimensions to be 

r = c ^!, c = 3 « 0.00171. (27) 

In Fig. 3 the classical spectral function and fits to a Breit-Wigner function are 
shown for various temperatures. From the fits one may extract estimates for the 
classical plasmon mass and width. This provides a possibility to compare non- 
perturbatively determined classical plasmon masses and widths with perturba- 
tive calculations and address the applicability of one-loop resummed perturbation 
theory at finite temperature. 

The calculations in the quantum theory can easily be carried over to the 
classical approximation used in the numerical calculation fl3 |. The one- loop gap 
equation for the classical mass parameter M c \ reads 



2 , a r d 2 P t 
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Figure 4: Mass M c \ (circles) and width T c i (squares) of the zero-momentum 
spectral function as a function of the temperature. The width is multiplied with 
25 for clarity. The data points are obtained from fits to a Breit-Wigner function, 
with the statistical error estimated from a jackknife analysis. The lines are the 
predictions for M c \ from the one-loop gap equation (full) and T c i from the on-shell 
two- loop contribution (dashed). 

The integral is logarithmically divergent in the ultraviolet, which reflects the 
logT contribution in Eq. (|25|) in the quantum theory. In the classical theory 
the divergence is regulated by the lattice cutoff. If desired, it can be matched 
(renormalized) by adjusting the mass parameter m in the effective classical theory 
[P~2]|. In order to compare the one-loop and the nonperturbative calculation we 
write the gap equation on the lattice, 

yj> ^ , 1 2 2 

M * = m 2 + —j- 2 £ ^— r^, p 2 = £^(l-cosap,), (29) 
A L ni,m P + IV1 cl i=i a 

with pi = 27Trii/Li, —Nij '2 + 1 < rii < Ni/2 and Li = aNi = L (i = 1, 2), and solve 
it numerically^ The result is presented in Fig. 4 for various temperatures. Maybe 
surprisingly we see that the nonperturbative determination of the plasmon mass 
and the resummed one- loop result differ at most a few percent, indicating the 
relative unimportance of higher-loop contributions for this quantity. Note that 
for our parameters the effective mass M c \ is small in lattice units (aM c \ ~ 0.3). 

An analytical expression for the width of the spectral function in the classical 
approximation can be obtained directly from the calculation carried out above. 

2 In the infinite volume limit the lattice gap equation can be written as = m 2 + 
\T/(4Tr)kF(n/2, k), with l/k = l + ^M^/i and F(ir/2, k) the complete elliptic function of the 
first kind. 



10 



As indicated, the perturbative width of the spectral function in the quantum 
theory is dominated by soft momenta and coincides therefore with the classical 
result to leading order: T c i is given by Eq. (p7[) after the replacement M — > M c \. A 
comparison between the perturbatively and nonperturbatively determined widths 
is presented in Fig. 4 as well. Again agreement between the two calculations can 
be seen, indicating that the dominant contribution is produced by the lowest order 
two- loop result. We emphasize that the width is a typical real-time quantity and 
therefore not easily accessible by other nonperturbative methods. 

5 Summary 

We have used the classical field approximation at high temperature and weak 
coupling to formulate a nonperturbative method for the calculation of spectral 
functions at finite temperature, with the help of the classical KMS condition. 
We focused on the one-particle spectral function in a scalar field theory in 2 + 1 
dimensions. For the temperatures investigated our numerical results indicate that 
the one-particle spectral function is a simple narrow peak: the spectral function is 
completely dominated by the plasmon. To compare with perturbation theory we 
calculated the one-loop resummed plasmon mass and the two-loop contribution 
to the width. Agreement between the nonperturbative numerical simulation and 
the perturbative expressions was found. These results provide a justification for 
the use of resummed perturbation theory for a weakly coupled scalar field in 
equilibrium as well as for kinetic approaches based on two-loop approximations 
close to equilibrium. 

It would be interesting to extend the analysis to more complicated spectral 
functions. Transport coefficients can be expressed in terms of spectral functions of 
composite operators. These are in general more sensitive to the ultraviolet scale 
and therefore to the Rayleigh- Jeans divergence in a classical limit. Neverthe- 
less, it might be interesting to test perturbative ideas against a nonperturbative 
numerical calculation. When sufficient care concerning the gauge symmetry is 
taken, methods applied here might also offer the possibility to gain further insight 
in nonperturbative aspects of hot gauge plasmas. 
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